Muon pair creation from positronium in a linearly polarized laser field 
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Positronium decay into a muon-antimuon pair by virtue of the interaction with a superintense laser 
field of linear polarization is considered. The minimum laser intensity required amounts to a few 
10 22 W/cm 2 in the near-infrared frequency range. Within the framework of laser-dressed quantum 
electrodynamics, the total reaction rate is calculated and related to the cross section for field-free 
electron-positron annihilation into muons. The muons are created with ultrarelativistic energies 
and emitted under narrow angles along the laser propagation direction. The dynamical properties 
of the muons are interpreted in terms of a classical simple man's model for the production process. 
We show that the most promising setup for an experimental investigation of the process employs 
two counterpropagating laser beams impinging on a positronium target, where the advantage of the 
coherent electron-positron collisions becomes evident. 

PACS numbers: 13.66.De, 36.10.Dr, 41.75.Jv 



I. INTRODUCTION 

Inside the most powerful laser fields available today [1] 
electrons and positrons can acquire kinetic energies in the 
GeV range. Next-generation laser devices are envisaged 
to generate even higher electromagnetic field intensities 
of ~ 10 23 W/cm 2 and above [2]. The GeV energy domain 
has also been demonstrated accessible via laser-wakcficld 
acceleration of electrons [3] and positrons [4] . In princi- 
ple, such high energies can be exploited to induce elemen- 
tary particle reactions like heavy lepton-pair or hadron 
production in e + e~ collisions. The energy scale associ- 
ated with superintense laser fields is thus more typical for 
particle physics than for atomic and molecular physics, 
where laser fields are employed traditionally. Apart from 
the high energies achievable, lasers can be utilized to 
generate well-controlled particle collisions at microscop- 
ically small impact parameters, which can lead to high 
luminosities [5-7] . Ions (protons) can also be accelerated 
by lasers [8], which may be used for triggering hadronic 
processes, such as pion production [9]. Consequently, 
intense laser fields offer alternative and supplementary 
ways towards high-energy physics [7, 10-13]. Similar ef- 
forts to merge laser physics with nuclear physics are be- 
ing undertaken successfully [14]. Laser-induced nuclear 
fusion [15], photofission [16], and neutron production [17] 
have already been observed in experiment, whereas other 
processes like coherent nuclear excitation [18] or muon- 
catalyzed fusion with laser assistance [19] have been dis- 
cussed on theoretical grounds. 

One of the most elementary processes in particle 
physics is the creation of a muon-antimuon pair from 
electron-positron annihilation, which has proven funda- 
mental for the understanding of other e + e~ reactions 
[20]. The threshold particle energy in the center-of-mass 
(cm.) frame required for this reaction is determined by 
the muon rest mass M 100MeV/c 2 and amounts to 
£" cm = 2Mc 2 . In conventional accelerator facilities the 
collision energy is provided by quasistatic electric and 
guiding magnetic fields. 



Motivated by the sustained progress in laser technol- 
ogy, we have addressed the question whether the process 
e + e~ —t is feasible with low-energy electrons and 

positrons in the presence of an intense laser field [12, 13]. 
The electrons and positrons are assumed to form initially 
a nonrelativistic e + e _ plasma or a gas of positronium 
(Ps) atoms [21]. Since the initial energy of the particles 
is far below 2Mc 2 , their annihilation into a muon pair 
cannot happen without the influence of the external field. 
We stress, however, that the muon pair is not produced 
by the laser field itself, as it has been achieved before with 
regard to e + e~ pair creation in the collision of a laser 
and an electron beam [22, 23]. In our case, this would 
require unrealistically strong laser fields of the order of 
the critical field = M 2 c 3 /eh for the muon. with the 
elementary charge e and Planck's constant h, which cor- 
responds to laser intensities / ~ 10 38 W/cm 2 . Instead, 
the muon pair is produced in an annihilating electron- 
positron collision, where the role of the laser field is to 
supply the required energy to the colliding particles. 

The case with initial Ps atoms is particularly inter- 
esting since it allows for coherent e + e~ collisions, which 
are characterized by microscopically small impact pa- 
rameters [7]. This is due to the fact that the electron 
and positron are initially confined to the atomic size of 
the bound Ps state and then, after instantaneous ioniza- 
tion, are coherently accelerated by the strong laser field, 
which drives the particles into opposite directions along 
the laser electric- field component. This leads to periodic 
e + e~ (rc)collisions whose energy is proportional to the 
laser intensity [24, 25]. Note that the magnetic field- 
induced ponderomotive drift motion into the forward di- 
rection is identical for the electron and positron due to 
the equal magnitude of their charge-to-mass ratios. In 
this respect the Ps system is distinguished from ordinary 
atoms where the heavy nucleus (or ionic core) stays be- 
hind, so that laser-driven collisions are suppressed at high 
field intensity [10]. In the case of strongly laser-driven Ps, 
the collisional impact parameter can be as small as the 
initial spatial separation of the particles which is deter- 
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mined by the Ps Bohr radius a>o w 1 A. This way, high 
particle current densities and collision luminosities are 
attainable [7, 13]. The combination of Ps atoms with 
intense laser fields might therefore be considered as an 
"e + e~ micro-collider". The threshold laser intensity to 
render muon production with significant probability from 
Ps or a low-energy e + e~ plasma is determined by the re- 
lation e\A\ > Mc 2 , with the laser vector potential A, be- 
cause the collision energy comes from the transversal mo- 
tion. For a near-infrared laser wavelength of A = 1 /im, 
the threshold intensity amounts to I = 5.5 x 10 22 W/cm 2 , 
which is almost reached by the most powerful present 
laser systems [1, 26]. 

The process e+e~ — > in a laser field has also 

been discussed within a two-step mechanism recently: In 
Ref. [27], an e+e - pair is assumed to be created in the 
collision of a laser pulse with a relativistic ion beam, 
and the muons are subsequently produced in a laser- 
driven e + e~ collision. The laser field strength in the 
ionic rest frame needs to approach the Schwinger value 
E$ = m 2 c 3 /eh = 1.3 x 10 16 V/cm (I ~ 10 29 W/cm 2 ). 
In this way, the well-established analogy between atomic 
ionization and e + e~ pair creation in strong laser fields 
has been extended to include also the recollision step. 
In a broader sense, the investigation of lepton-lepton in- 
teractions in laser fields has a long history, comprising 
laser-assisted e"e" (Tvfoller) [28], e+e~ (Bhabha) [29], 
and e~ \i~ scattering [30]. In contrast to the process 
under consideration here, these lepton-lepton scattering 
processes can also take place without a background laser 
field, and the presence of the latter merely modifies the 
field- free cross-section. 

In the present paper, we give a comprehensive account 
of the reaction Ps — ► n + in a laser wave of linear polar- 
ization, this way extending our recent study in Ref. [13]. 
The Ps atom is supposed to be initially at rest and in 
its ground state. The pair is created in the laser- 

driven e + e~ collisions as described above (see also Fig. 1). 
On the one hand, we perform a quantum electrodynam- 
ical calculation to obtain the total production rate and 
muon spectra; both analytical and numerical results are 
presented. The occurence of sums of generalized Bessel 
functions of very high order > 10 10 constitutes a major 
difficulty of the case of a linearly polarized driving field. 
On the other hand, particular emphasis is laid upon the 
derivation of a semi-classical model for the muon pro- 
duction process, representing one of the main extensions 
of [13]. The model allows to interpret the typical muon 
momenta and the total process probability in a simple 
and intuitive manner. It is proven by a saddle-point con- 
sideration and resembles the famous three-step model of 
laser-driven recollisions [25] , but also exhibits character- 
istic differences. The saddle-point treatment moreover 
provides us with a suitable asymptotic expansion to han- 
dle the high-order generalized Bessel functions. Further- 
more, a comparison is drawn with muon pair creation 
from Ps in a laser field of circular polarization, which was 
treated in Ref. [12]. In this situation it was found that 



the muon production is essentially suppressed, because 
the e + e~ rescattering occurs at large impact parameters. 
The latter are determined mainly by the classical electron 
and positron trajectories in the circularly polarized field 
and have the macroscopic dimension g <~ A£ <~ 10~ 2 cm, 
where £ ~ 200 is the laser intensity parameter [see Eq. (6) 
below]. Instead, in a linearly polarized field the colli- 
sional impact parameter can be microscopically small, 
g ~ ao w 10~ 8 cm, which enhances the probability for Ps 
decay into a muon pair by several orders of magnitude. 
Finally, the experimental feasibilty of laser-induced muon 
production is discussed in detail. It is shown that an ob- 
servation of the process will be rendered possible in a 
crossed-beam setup employing high-power laser sources 
of the next generation combined with high-density Ps 
samples. 

The paper is organized as follows. In Sec. II we de- 
rive general expressions for the probability of laser-driven 
Ps decay into muons in two ways. First, we proceed in 
Sec. II. A along the lines of our earlier paper [12] and 
evaluate the process amplitude by virtue of a Fourier 
expansion, where the Fourier coefficients are given by 
generalized Bessel functions. In Sec. II. B we employ an 
alternative approach based on the saddle-point integra- 
tion method, which will prove useful in various respects: 
We show that important features of the considered muon 
production process can be understood in classical terms 
(see also Sec. III.B) and derive an asymptotic expansion 
of the generalized Bessel function [see Eq. (30)], which 
later serves us in our numerical computations. These spe- 
cial functions frequently occur in calculations of intense 
laser- matter interactions [31-33]. Based on the general 
expressions at hand, the relevant physical features of the 
process are elaborated in Sec. III. The kinematics of the 
produced muons is analyzed in Sec. III. A and explained 
in Sec. III.B within an intuitive semi-classical model. In 
Sec. III. C we derive a closed analytical approximation to 
the total creation rate [see Eq. (54)], which exhibits the 
main physical features of the process and can be inter- 
preted in terms of the cross section for the field-free re- 
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FIG. 1: Feynman graph for muon pair creation from electron- 
positron annihilation in a background laser field. The ar- 
rows are labeled by the particle's free momenta (p±, P±) out- 
side and the effective momenta (q±, Q±) inside the laser field. 
The virtual photon has four-momentum q. The electron and 
positron are assumed to form a Ps atom initially. 
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action e + e 



and the quantum mechanical dis- describes the coherent average over the amplitude 



pcrsion of the initially bound e + e~ wave packet in the 
laser field. Afterwards, we briefly consider in Sec. III.D 
the related process of muon pair production by a su- 
perstrong laser wave of linear polarization which inter- 
acts with a nonrclativistic e + e~ plasma. The influence 
of the laser polarization is revealed in Sec. III.E. In the 
following Sec. IV we present results of numerical calcula- 
tions, compare them with previous ones for circular laser 
polarization, and discuss their relevance for an experi- 
mental investigation of laser-driven muon creation. Here, 
also more complex driving field configurations and their 
potential advantages are described. We summarize our 
main results and conclusions in Sec. V. 

Relativistic units with K = c = 1 are used through- 
out, except where otherwise stated. The finestructure 
constant is a — e 2 w 1/137. We employ the metric 

tensor g^ u = diag(H ), so that the scalar product 

of two four- vectors p M = {p°,p) and g M = (q°,q) reads 
(pq) = p°q° — pq. Furthermore, Feynman slash notation 
i> = ilP) is used, with the Dirac matrices 7 P . 



II. THEORETICAL FRAMEWORK 

A. Transition amplitude and reaction rate 

The rate for positronium decay into muons in a 
strong laser field can be calculated within the frame- 
work of laser-dressed quantum electrodynamics (QED). 
The driving laser field is assumed to be a monochromatic 
plane wave of linear polarization, described by the clas- 
sical four-potential 



A^(x) = e^acos(kx), 



(1) 



with the polarization four-vector e M = (0,1,0,0), the 
wave four-vector fc M = u)(l, 0,0,1), and the amplitude 
a. In the following consideration, the laser parameters 
have typical values of ui ~ 1 eV and a > 100 MV, cor- 
responding to a high intensity I > 10 22 W/cm 2 in the 
near-infrared frequency domain. When an electron and 
a positron are exposed to such an intense laser wave, 
the interparticle Coulomb interaction is much weaker 
than the coupling to the external field and may be ne- 
glected. Within the spirit of the strong-field approxima- 
tion [31, 34], the amplitude of the process Ps — > 
can therefore be expressed as a superposition integral [12] 



>Ps 



Vvl 



' rp *(p)5 e+e 
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(P), 



(2) 



with a normalization volume V and the Fourier transform 
$(p) = 8y/Tra 3 /2 1 (1 + alp 2 ) 2 of the Ps ground state of 
Bohr radius ao- The latter function is also known as 
the atomic Compton profile. The Ps atom is initially 
assumed at rest, and p denotes the relative momentum 
of the electron-positron two-body system. Equation (2) 



Se+e-(p) = -ia / d 4 x J d 4 y* p+ , s+ (x)^ p _ !S _(x) 
xD„ v (x - yj*p_,s-(yh v *p + ,s + (y) (3) 

for the process e + e~ — ► in a laser wave (cf. Fig. 1). 

The initial electron and positron momenta are p± = ±p. 
Their wave functions in Eq. (3) are given by the Volkov 
states [35, 36] 



(4) 



with the phase 

(kp±j 



e 2 a 2 
8(kp±) 



sin[2(fcr)]. 



In Eq. (4), p± = (po,±p) arc the initial free four- 
momenta of the electron and positron (outside the laser 
field), s± denote the particle spin states, the u p± , s± are 
free Dirac spinors [37], and 



A{kp ± ) 



(5) 



are the effective four-momenta of the particles inside the 
laser field. In accordance with Eq. (5), the momentum 
vectors satisfy the relations q± = p± and q+ + q^ = 0, 
where the label _L denotes the component perpendicular 
to the laser propagation direction. The corresponding 
effective mass reads to 2 = q\ = m 2 (l + £ 2 ), with the 
dimcnsionless laser intensity parameter 



(6) 



The Volkov states in Eq. (4) are normalized to a num- 
ber density of one particle per unit volume, and to a 5- 
function in q± space [36, 38]. Analogous expressions hold 
for the Volkov states v Pp ±i s±j the free momenta P^, the 
spin states S± , the effective momenta Q± , and the effec- 
tive mass M* = (M 2 + m 2 £ 2 ) 1 ^ 2 of the muons. Moreover, 
we apply in Eq. (3) the free photon propagator 



D^(x - y) 



(2^) 4 q 2 9 ' 



(7) 



We comment on various approximations that are im- 
plicitly contained in the above equations. (1) The ini- 
tial spin coupling in the Ps ground-state is neglected in 
Eq. (2). It is of minor importance in the presence of the 
strong laser wave; the corresponding hyperfine splitting 
between the ortho- and para-states amounts to less than 
1 meV which is negligibly small on the energy scale set 
by the external field. Moreover, later on we will aver- 
age over the spin states [see Eq. (22)] so that our ap- 
proach applies to a statistical ensemble of Ps atoms. A 
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spin state-resolved calculation could be performed by in- 
serting appropriate projection operators [39]. (2) The 
amplitude in Eq. (3) fully accounts for the interaction of 
the leptons with the laser field, while their interaction 
with the QED vacuum is taken into account to lowest 
order, in accordance with the Feynman graph in Fig. 1. 
The interplay between the Coulomb field and the elec- 
tromagnetic wave during the very beginning of the laser 
field's rising edge might lead to small distortions in the 
momentum distribution (3) Our treatment more- 

over neglects the radiative energy loss of the laser-driven 
electron and positron due to Thomson scattering, which 
is rather small for the laser parameters under considera- 
tion [12]. (4) Finally, we ignore in Eq. (7) modifications 
of the photon propagator due to the presence of the ex- 
ternal field. These modifications become important only 
at substantially higher laser intensities / > 10 29 W/cm 2 
[40, 41]. 

The space-time integrations in Eq. (3) can be per- 
formed by the standard method of Fourier series expan- 



sion, using the generating function of the generalized 
Bessel functions [31]. The latter can be expressed via 
ordinary Bessel functions according to [31] 



J„(U,V)= Jn-2i(u)M v )- 



(8) 



For the amplitude, we obtain this way 
Se+e-(P) = -i(27r) 4 a 



^Q° + Q°- 

j y ~ at 



n,N 



x 5(q+ +q--q-nk) S(Q + + Q_-q- Nk), 

(9) 

with the spinor-matrix products 



(i 

\2{kp+) 

M,(P + ,P.\N) = u P _,s.[l»B N (gfe - 



e 2 a 2 k^ft 



2(kp-)J 2{kp+){kp-) 
earuM \ e 2 a 2 k^ , 

2{kPT) ) ° n ~ 2{kP + ){kP.) Dn < UP +> S +- (10) 



r 



The coefficients of the electronic spinor-matrix products 
read 

b n = J n {u,v) , 

Cn = ^ [J n -l(u, V) + J n+1 (u,v)] , 

d n = 2[J n -2(u,v)+2J n (u,v) + J n+2 (u,v)} 

with the arguments 

ea(ep-) ea(tp + ) 
(hp-) (kp+) 



e 2 a 2 r 



1 



(hp-) (kp+) 



(11) 



Note that v ~ v = — to£ 2 /(2w) is practically constant 
for the relevant values of |p| ~ ma. The coefficients B^, 
Cn and Dn of the muonic spinor-matrix products have 
the same structure, with the arguments given by 



U 
V 



ea(eP_) ea(eP+) 
~JkP~7) (WVf' 



e 2 a 2 



1 1 

+ 



(fcP_) (kP + ) 



(12) 



According to the energy-momentum conserving 5- 
function at the first vertex of the Feynman graph, the 



integer number n in Eq. (9) counts the laser photons that 
are emitted (if n > 0) or absorbed (if n < 0) by the elec- 
tron and positron. Similarly, N is the number of photons 
exchanged with the laser field by the muons. Denoting 
the total number of absorbed laser photons by r := N — n 
and integrating over the virtual photon momentum yields 



Se+e-(P) = -i(27T) 4 a 



M 



xY,M^(p+,p-\n)M,{P + ,P-\n + r) 



,%4 



Q-+rk) 



(q+ +q-- nk) 2 



(13) 



By Eq. (8) and the properties of the ordinary Bessel func- 
tions [42] , the main contribution to the amplitude comes 
from photon numbers n rts 2v rts 2v = —m{; 2 /u). Conse- 
quently, the square of the virtual photon momentum 

(q+ + q--nk) 2 = 2[m 2 + (q+q-) - n{kp+) - n(kp-)\ 
w 4(m 2 - m) w 8m 2 £ 2 (14) 

is practically constant and may be pulled out of the sum- 
mation over n in Eq. (13): 

M 1 



S e +e-{p) 



-i(27r)V 



rn 



q° + q- 



Q+Q 



o no 8to 2 £ 2 



■5 



x Y,^(p + ,p-\n)M^P+, P-\n + r) 

n,r 

x6(q+ + q_-Q+-Q_+rk). (15) 

A particular implication of Eq. (14) is that the process 
proceeds nonresonantly (i.e. q 2 ^ 0) [43]. For later use 
we note further that the typical energy of the virtual 
photon is 



(jo = <7+ + 9_ — nuj w 2m£ 2 . 



(16) 



The sum over n in Eq. (15) can be performed analyti- 
cally by virtue of an addition theorem for the generalized 
Bessel functions [31] 

oo 

J2 J n (u,v)J n+r (U,V) = J r (U-u,V-v). (17) 

n— — oo 

In this manner we find 

5^M"(p + ,p_|n)M M (P + ,P_|n + r) = J r u^ 

n 

+J£ r (u'% + + L r (u"W„ + v*% + w»U^) 

+ Mr^W^ + 10*%) + NrW^W^. (18) 

Here we have introduced the notation 



— ^lp+,s+'~f^ L ^p— .S — 



w M = u 



e 2 a 2 fc^ 
'2(fcp+)(fcp_) 



and similarly U^, V^, and W M for the muons. The coeffi- 
cients in Eq. (18) read 

J r = J r (U — u,V — v), 

K r = - (J r _i + J r +i) , 

L r = — (J r -2 + 2J r + J r +l) , 

M r = - ( J r _3 + 3 J r _i + 3 J r+ i + J r +3) , 

5 

Af r = — ( J r _4 + 4J r _ 2 + 6J r + 4J r+2 + Jr+4) • 

After these transformations of the amplitude for the 
process e + e~ — ► in Eq. (3) we come back to the 

amplitude for the decay Ps — ► in Eq. (2). We need 

to multiply Eq. (15) by the Compton profile $(p) of the 
Ps ground state and integrate over the relative momen- 
tum p. This complicated integral can be performed in 
an approximate way as follows: 

Within the momentum range given by <fr(p), the elec- 
tronic spinor-matrix products in Eq. (19) may be con- 
sidered as constant and pulled out of the integral. The 
same holds for the kinematic factors q± w to£ 2 /2 and the 
energy-momentum conserving 5-function. (The latter is, 



in fact, a subtle issue and dealt with in Appendix A). 
Hence, we are left with integrals of the form 



u, V — v). 



(19) 



Despite the small p values contained in the Compton pro- 
file, the high-order Bcsscl function [cf. Eq. (36) below] in 
the integrand is heavily oscillating, which leads to a small 
value of J r . The oscillatory damping of the amplitude is 
due to a destructive interference of the various partial 
waves within the Ps wave packet, which can equivalently 
be described as wave packet spreading [see Eq. (57)]. The 
evaluation of the integral (19) is performed in Appendix 
B. With AV = V-v , the result reads 



2 9 



2/3 



(a£AV) 



1/3 



-J r (U,AV). (20) 



/ 7t^Gq \ma 2 t; J a 2 |t>o| 
Now we can write the square of the Ps decay amplitude 



as 



|5 Ps | 2 = (27r) 4 a 2 



m 2 M 



1 



q^q - Q° + Q°- 64m 4 £ 4 



E 



J2M^ P+ ,p-\n)M^P + ,P-\n + r) 



X S(q+ + q_ - Q+ - Q_ + rk)VT , 



(21) 



where the symbol ^ n indicates both the sum over n and 
the average over the Ps ground state, as described above. 
p± and q± in Eq. (21) are to be understood as some typ- 
ical values of the electron and positron momenta. The 
factors of volume V and time T come, as usual, from 
the square of the 5-function. Note that the energy- 
momentum conserving 5-function implies + Q^z = 0. 
From Eq. (21) we obtain the total rate for laser-driven Ps 
decay into muons by averaging over the initial spin states, 
summing over the final spin states, and integrating over 
the final momenta: 



R 



Ps 



1 [<PQ + fd 3 Q_ 1 ^ 



s±,S± 



Based on Eq. (22), we have performed numerical calcu- 
lations; the results will be shown in Sec. IV. It is also 
possible to evaluate the integrals in Eq. (22) in an ap- 
proximate way by analytical means. This derivation is 
performed in Sec. III.C. 



B. Saddle-point approach 

Due to the periodicity of the laser field, the exponen- 
tials in the amplitude (3) are highly oscillating functions. 
As an alternative to the series expansion into general- 
ized Bessel functions as performed in Eq. (9) , one can 
evaluate the amplitude by applying the saddle-point in- 
tegration method. Our aim here is twofold: on the one 
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hand, to derive a simplified expression for the Ps decay 
amplitude (21), and on the other hand, to interpret the 
saddle-point integration results via a simple-man model 
which will lead to an intuitive interpretation of the pro- 
cess probability (see Sec. III. B). 

It is convenient to introduce light-like coordinates 
x± = x° ± x 3 and y± = y° ± y 3 , so that the inte- 
grals over dx + , dx- 1 , dy + and dy 1 - all result in <5 func- ^as ^ e solution 
tions. Among them are d(q+ + q^z — q^) = S(q ) and 



where AQ = Q Q + + Q_ — q . The saddle point equation 



-p- = AQ - ULucos(ujy-) - 2Vu cos(2wy_) = (27) 

OA)— 



%+ -q 3 + + q°--q 3 --q° + q^) = S(2p Q - q° + q 3 ), which 
we exploit to integrate by dq ± and dq 3 . This yields 



= — la-K 



dx- 



dq 



Po(qo~Po) 



— i(q^_-\-q^_ —qo)x^-\-iu s\n(iox_ )-\-iv s'm{2ux _ ) 

AC ' 

x gi(Q++Q_ -Qo)y- -iU sin(w2/_ )-iV sin(2u>2/_ ) 

x^(p + ,p_;i_)^(P +I F_;^ 
x5(Qjr + Q_:)S(Q° + -Q 3 + + Q°_- Q 3 _ - 2 Po ) 

(23) 



with 



M^(p+,p-;x-) 



m 



q° + q°- 



2{kp+) 



r 



i - 



2(kp- 



.(24) 



and accordingly _M M (P + , P_; y_) for the muons. The 
integrals over <ix_ and dy_ can separately be done 
by the saddle-point method. We explicitly demon- 
strate the y_ integration by considering first the 
term in .A/f M (P+, P_; y_) which is proportional to 
wp_,s_7'*Mp + ,s + and independent of y_. The corre- 
sponding integral is denoted as 



/+oo 
-OO 



dy- e 



if(v-) 



(25) 



with the phase given by 

= AQy_ - E/sin( W </_) - Vsm(2ujy-), (26) 



U 



cos(wy_) = -gp7± 



_CP_ 

64V 2 



+ 



AQ 1 

4Vlo 2 ' 



(28) 



The resulting saddle points y_ indicate the regions giving 
the main contribution to the process. Physically speak- 
ing, they determine the laser phases when the muons are 
created predominatly. The corresponding saddle points 
of the x- integral represent, accordingly, the laser phases 
when the e+e~ annihilation occurs. This fact will be uti- 
lized in Sec. III.B below. 

Due to the periodicity of the phase (26), we can rewrite 
Eq. (25) as 



au e i(U sin 9+V sin 20-U0) 



' 277 N 



T 1 ^ r-Z 

1 +°° f-2TT 

±_ c -27riW^ / jQ/ c i(U sin6' + V sm2d'-vd') 

+ 0O 

2tt S ( A Q - Nuj)J n (U,V) , (29) 



N=-co 



with v ee AQ/lj, 9 ee -lux- and 9' = 9 - 2ttN. In 
the last step, the integral representation Jn(U, V) = 
(1/2tt) J Q 27r exp[i(C/ sin 9 + V sin 29 - N9)]d9 from Ref. [31] 
has been applied. If the 9' integral is evaluated instead 
with the help of the saddle points in Eq. (28), one obtains 
the asymptotic expansion 



1 /cos(iVj ?+ -C/sin?7 + -ysin27 7+ + a(+)f) cos(7Vr ? _-C/sinr ? _-ysin2?7_-a(-)f) \ 

Mu ' v) = ^p { 7W^T\ + 7W^r\ ) m 



provided that the saddle points arc real and isolated from _ j 

T}-$- — arccos i ~ __t p 

each other. Here we have introduced the notation V 8V 



U_ 

8V 

CT (±) _ sgn(y sin?7±). 



U 2 N 1 The asymptotic representation (30) has been used in 

04y2 ' 4T/ 2 ' our numer i ca l calculations (see Sec. IV) . It considerably 
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simplifies the evaluation of high-order generalized Bessel 
functions in comparison with the series expansion (8). 
For the reader's convenience, the details of the deriva- 
tion are given in Appendix C. The formula (30) is also 
contained as Eq. (3.8) in the comprehensive study of 
asymptotic expansions of generalized Bessel functions in 
Ref. [32]. Its application presumes real parameters U, V, 
and N of large magnitude. The range of applicability is 
specified more precisely in the appendix. 

We briefly check that the present approach leads to the 
same result as in Sec. II. A, where first the integrations 
over d 4 x and d 4 y have been performed and afterwards 
the d 4 q integral. By using Eq. (29), the expression 



, 4 / d 10 

S e +e- ~ -lain / — - - 

J po{qo -Po) +ie 



S(Aq — nuj) 



xJ n (u, v) J2 S ( A Q- Nlu)J n (U, V) 

N 

m M 



x 



<#g° JQ° + Q° 



xS(Q°+ -Q 3 + + Q°_ -Q-- 2po), (31) 

with Aq = q^_ + q°_ — qo , is found for the leading term in 
Eq. (23) which contains the integral (25). The integration 
over dqo in Eq. (31) can be performed with the help of 
the first S function yielding 



n,N 



J n (u,v)J N (U,V) m M 
(q + + g_ - nk f y^JT ^/ Q o go 

x 5 {A) (Q+ + Q_ - q+ - q- - (N - n)k) . 

(32) 



Here we used the relation 



Po(qo-Po) ~ m(q1 + q°_ - nw - m) 



nkf 



which follows from Eq. (14), and combined the remaining 
5 functions in Eq. (31) into a four-dimensional one by 
virtue of the identity 



2m 



(N- 
q+,z + q- 



n)to 

z + (N-n)tu. 



We note that the structure of Eq. (32) coincides with that 
ofEq. (13). 

So far, we have only dealt with the space-time inde- 
pendent parts of 7W i (p + ,p_; xJ) and M. fl {P + ,P-\yJ) 
in Eq. (23). There are also terms containing cos(wy_) 
and cos 2 (wy_) in A4^(P+, P_; y_) which can be handled 
in a similar manner as before. For example, by writing 
cos(wy_) = i(e lwy - +e _lwy -), the integral corresponding 
to Eq. (25) becomes 



dy- 



if-iv-V 



with the phases 

f±(y-) = (Q° + + 0°_- ( z ±^)y--[/sin(^_) 
+2Vsin{2ujy-). 

By recapitulating the previous calculation, we see that 
the new phases f± lead to the replacement of AQ by 
AQ ± u) or, correspondingly, of N by N ± 1. Thus, we 
obtain 

+ oo . 

l 1 = 2n S(AQ-Nlo)-[J n+1 (U,V) + J n ^(U,V)] 

N=—oo 

in analogy with the Bessel function treatment in Sec. II. A 
[cf. Eq. (11)]. As expected, the saddle-point approach is 
equivalent to the treatment which employs a Fourier ex- 
pansion into Bessel functions, when a suitable asymptotic 
representation of the latter is applied. For the analogous 
result in the case of a circularly polarized laser field, we 
refer to Ref. [44] . 



III. ANALYTICAL DEVELOPMENTS AND 
RESULTS 

A. Kinematical analysis 

By combining energy-momentum conservation as ex- 
pressed by the 5 function in Eq. (21) with Bessel function 
properties, we can gain rather comprehensive information 
on the kinematics of the produced muons (see [12] for a 
similar analysis). The energy-momentum conservation 
condition 

4m 2 + 4vu)ui w (q+ + q- + rk) 2 = (Q+ + Q-f> 4M 2 

implies that the minimal number of laser photons to be 
absorbed is 



M 2 - m 2 _ M 2 
turn torn 



(33) 



It is interesting to note that even this minimal pho- 
ton number gives rise to ultrarelativistic muon momenta 
Q± w q± +r m j n a>/2 3> M*. The 5 function further yields 



Q°+ - Q+,z +Q--Q- 



q+ - q+,z + q_-q- 

P° + -P+,z+P- -P- 



= 2p a w 2m , 

which implies 

(kP+) w (fcP_) w (kp+) w (kp-) w cum (34) 

under the symmetry assumption Q°_. The latter 

condition will always be assumed in what follows, which 
allows us to simplify the notation by dropping the parti- 
cle indices ± from time to time. From Eq. (34) we obtain 
in particular the relation 



V w v 



(35) 
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between the Bessel function arguments of Eqs. (11) and 
(12). By expressing the virtual photon momentum in 
terms of the muon momenta 



8m 2 £ 2 



(q+ + q- - nkf = (Q+ 
4M 2 + AQl - ANujm , 



Q_-Nkf 



we find that the typical number r — N — n of absorbed 
laser photons is 



M 2 + Q 2 ± 



(36) 



Due to Eq. (35), the difference AV = V — vq has a smaller 
magnitude than V or vq themselves and, consequently, 
than r. In order to obtain non-ncgligibly small values of 
the Bessel function J r (U, AV) in Eq. (20), it is therefore 
neccessary to have U w r. By Eq. (36) this implies 



P x = V2m£ (l + v 7 ! - k 2 ) 



(37) 



with k 2 = (M 2 + P 2 )/(2m 2 £_ 2 ). Requiring k 2 < 1, we 
find 



\P y \<P^ = V2my/e-^, 
along with the minimal laser intensity parameter 
c. = M 

smin — 



V2' 



(38) 



(39) 



We emphasize that £ m i n is minimal in the sense that it 
is the smallest £-value which allows for a non-negligibly 
small process pobability. In a background laser field the 
decay Ps — > is energetically possible at any laser 

intensity in principle, since the external field represents 
a practically infinite energy reservoir. However, only for 
laser intensities satisfying Eq. (39) the probability of this 
reaction can be significant, as it guarantees a sufficiently 
efficient photon absorption. Note moreover that Eq. (39) 
confirms the naive estimate of the threshold intensity 
given in the introduction. This is, in fact, remarkable 
since the muons have to be created in the field with their 
laser-dressed mass M*, which is significantly larger than 
their bare mass M (see also Sec. IV. B). 
By combining Eqs. (36)-(38) we find 



1 . M 2 + 3m 2 £ 2 
Q z « 2 + + ™> M 2^ 



which translates into 



2 £2 



M 2 + 2m 2 ^ 
2nT 



(40) 



(41) 



by Eq. (5). Close to threshold (£ ss £ m i n ), the produced 
muons have typical momenta of 



M, Py^O, P Z 



M 2 
m 



(42) 



The total number of absorbed laser photons is r w 
2M 2 /uim, which is of the order of 10 10 at u> = 1 eV and in 
agreement with the energy conservation law outside the 
laser field: rui ss 2(P — m). Note that, even at threshold, 
the produced muons are highly relativistic such that their 
life time in the lab frame is considerably increased due 
to time dilation. The typical muon momenta in Eq. (42) 
can be interpreted in terms of a classical simple man's 
model as described in the following Sec. III.B. 

For later use we note that the typical muon momenta 
given above imply, that the particles are emitted into a 
narrow cone along the laser propagation direction. In- 
deed one can show [12], that the opening angle of this 
cone is 



2m 



(43) 



The occurence of a maximum polar emission angle re- 
sults from the fact that the electron and positron move at 
relativistic longitudinal velocity in the laboratory frame. 
In the cm. frame, the muons are emitted isotropically. 
In contrast to the rotationally symmetric case of circu- 
lar laser polarization [12], there is also a maximum az- 
imuthal angle. According to Eqs. (37) and (38), we have 
P™ ax < P x and thus P x k, P ± . Consequently, 



arcsin 



pmax 
V 

Pi 



, pmax 

y 

P± 

< 


(for 


pmax 
V 


7T 

' 2 


(for 


pmax 
V 


pmax 


/- 


e 2 - 




e 



« p±), 



p±) 



(44) 



The angle </> max determines the maximally allowed devi- 
ation from the laser polarization axis. 



B. Simple man's model of muon creation 

By using a classical simple-man's model, one can give 
an intuitive interpretation of the typical momenta in 
Eq. (42). This is similar to the famous three-step model 
of above-threshold ionization and high-harmonic gener- 
ation in strong laser fields, which can also be derived 
within a saddle-point approach by equipping the saddle 
points with a classical meaning (see [45] for a review). 
In this simplified model the electron and positron arc 
released from the bound Ps state by a strong laser field, 
excurse in the laser field as free particles acquiring energy 
from photon absorption, and produce a muon pair at the 
moment of the recollision. Let us consider the process at 
threshold (£ = £ m i n ), where one may assume that the par- 
ticles are created with zero momentum in the (primed) 
electron-positron cm. frame, whenever the laser vector 
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potential reaches its maximum amplitude. We note that 
the latter implies vanishing electric field strength of the 
laser. This agrees with the intuitive classical picture of 
laser-driven e+e~ (rc)collisions, which always occur near 
zero electric field strength. It moreover indicates that the 
muon pairs are not produced directly by the laser field 
but rather by the energetic electron-positron collisions. 

First, we demonstrate that the above assumptions in- 
deed lead to Eq. (42). Let P '(to) = 0, where To denotes 
the laser phase when the production occurs. The solu- 
tion of the classical equations of motion with this initial 
condition reads 



= e[A(r)-A(r )}, 



P'Ar) = —[A(r)-A(r )f 



(45) 



where the laser vector potential is given by A(t) = 
a cost, with the Lorentz-invariant phase t = uu{t — z) 
[see Eq. (1)]. After having left the laser field, the momen- 
tum components of Eq. (45) become P' x = M, P' y = 0, 
P' z = M/2. Here, the conditions A(t ) = a and £ = £ min 
have been used. The Lorentz transformation to the lab 
frame now yields the momentum components of Eq. (42), 
with the reduced velocity of the frame transformation 
given by Eq. (48) below. 

The assumptions that our simple man's model is based 
upon, can be infered from the saddle-point treatment of 
Sec. II. B. With regard to the electron-positron annihila- 
tion vertex, the saddle points X- are determined by 



cos(wx_) = 



u / r 
± V Mv 2 



n 1 
Av 2 



(46) 



Since |u| < ^am/uo <C \v\ m £ 2 m/cj and n s=a 2v in our 
situation, Eq. (46) becomes cos(wx_) w ±1. This means 
that the annihilation occurs at the laser phases x_ when 
the amplitude of the vector potential (I) attains a max- 
imum. On the femtosecond time scale of the laser oscil- 
lation, the muon creation occurs directly afterwards, in 
agreement with one of our model assumptions. This cir- 
cumstance moreover implies that, at £ = £ m in, the muons 
are generated at rest in the cm. frame, as we show now. 
The classical energy of a particle of charge — e in a laser 
field given by Eq. (I) reads 



(eP) 

Po(y-) = Po - eac^y— y cos(wy_) 



2 2 

e a u> 
2(kP) 



cos 2 (wy_) . 



where P M is the initial momentum and y_ = t — z. Con- 
sequently, the saddle point equation (27) can be recast 
in the form 



P°( 2/ _)+P ( y _)= <Z o, 



(47) 



where qo ~ 2m£ 2 according to Eq. (16). We need to boost 
this equation into the cm. frame. At maximum value 



of the laser vector potential, the reduced velocity of the 
electron (or positron) along the z axis is 

Pz(x-) e 



i - 



(48) 



where po(x~ ) and p z (x- ) are the classical energy and mo- 
mentum of the electron. This corresponds to a Lorentz 
factor of 7 ~ £/ \/2 and yields 



7(9o - Pq z ) ~ 7 



1 



1 



{q - 2p ) 



r, |!L + 2m 7 sa 2V2mi , 
2 7 

where q z denotes the z-component of the virtual photon 
momentum . At the threshold intensity [cf. Eq. (39)] wc 
therefore obtain for the sum of the muon energies 

Pl(y-)' + P°{y-)' = q' ~ 2V2mCmin = 2M , (49) 
i.e., the muons are created at rest. 

C. Analytical estimate of the total rate 

From Eq. (22), one can derive by analytical means an 
approximate formula for the total rate of muon produc- 
tion from a single laser-driven Ps atom (see also [12]). 
To this end, we consider the contribution to the rate 
stemming from the first term on the right-hand side of 
Eq. (18): 

M 2 f d 3 Q+ [ d 3 Q- 



(2tt) 2 2 6 m 4 C 8 



X^\u»U^6(q + +q 

s±,S± 



Q+ J Q°- 

Q+-Q-+rk).{50) 



Note that, when squaring Eq. (18), eighty other terms 
of a similar structure arise. According to our numerical 
calculations, the term accounted for in Eq. (50) gives a 
major contribution to the total rate. With the help of 
the four-dimensional 6 function, we can integrate over 
d 3 Q+dQ°_ and find 



R 



Ps 



a 2 M 2 Ml 
(2tt) 2 2 8 m 6 e 8 

r s±,S± 



dcos 9q_ 



(51) 



From the Jacobian of the energy-conserving 8 function, a 
factor « [2(1 — cos^q_)] _1 ss (M*/2m) 2 arises here. We 
point out that this factor has been overlooked in [12], 
however, without affecting the main conclusions of this 
paper. The spin sum in Eq. (51) can be converted into a 
product of two traces in the usual way, 



8±,S± 



Tr I y ^-+ m 7 ^+- m 
2m 2m 



, P_ + M f+ 
xTr 7 / — 7, 



AT 



2M 



2M 
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The standard trace technology yields 
2 



%iU — 



[(p_P_)(p + P + ) + (p_P + )( P+ P_)] 



m 2 M 2 

, 2(P+P_) , 2(p+p 



M 2 



+ 



According to the kinematical analysis in Sec. III. A, a typ- 
ical value of T u u is 



16m 2 £ 



2f4 



uU ~ M2 , 



(52) 



which can be pulled out of the integration in Eq. (51). 
We proceed by performing the further approximations 



^[J r ((7,AT/)] 2 « 1, AV 1 / 3 ^ \v \ 1/3 , (53) 



Due to relativistic time dilation, the oscillation period in 
the cm. frame is largely enhanced so that the recollision 
time at £ 3> 1 equals i£ cc ~ 7T7 /lu, with the Lorentz factor 
7 <~ £ [see Eq. (48)]. Consequently, at the collision the 
wave packet has spread to a size Ax'Ay'Az' ~ (a£A) 3 . 
The corresponding particle current density leads to a re- 
action rate of Rp s <~ er/(a£A) 3 . By boosting this rate 
into the lab frame we obtain 



£«A)3 



(57) 



which agrees with Eq. (54) up to factors of order unity. 
This demonstrates that simple rate estimates for (nonres- 
onant) nuclear or particle reactions in strong laser fields, 
which are based on semi-classical arguments and field- 
free cross sections, may be considered reliable (see, e.g., 
[5-7]). 



where the maximum polar and azimuthal angles are given 
by Eqs. (43) and (44). The first two estimates in Eq. (53) 
mean that the production rate is assumed as approxi- 
mately constant within the mainly contributing angular 
domain. The factor of 4 in the 4>q_ integral is due to the 
fact that the integration range of interest extends over 
(-^max, 0max) and (it - (f) maK , 7T + <j) max ). The third re- 
lation would exactly be true [31], if the muon momenta 
and, thus, U and AT^ were independent of r. The last es- 
timate in Eq. (53) holds since AV = V — v and Vq have a 
similar order of magnitude; i.e., in general V and v w vo 
in Eq. (35) compensate only partially, as is corroborated 
by our numerical calculations. As a result, we arrive at 
the formula 



R 



a 



7T 2 m 2 £ 2 ao 



which is the desired rate estimate for Ps decay into muons 
induced by a linearly polarized laser field. The analyti- 
cal approximation (54) is confirmed by direct numerical 
evaluation of Eq. (22), as shown in Fig. 2 in Sec. IV. 

Equation (54) has an intuitive interpretation in terms 
of the cross section for the field-free process e + e~ — > 
[20] 



a = 



Att a 2 



'1- 



4M 2 

cm 



i + 



2M 2 

cm 



(55) 



and the electron-positron wave packet size, which is grow- 
ing in the laser field due to quantum mechanical dis- 
persion. In our case, the cm. energy in Eq. (55) is 
E cra fa 2\f2m^ according to Eq. (14). We note that 
both Eq. (54) and (55) vanish at the intensity (or energy) 
threshold. The wave-packet spreading can be estimated 
as follows. In the (primed) electron-positron cm. frame, 
the initial momentum spread is Ap' <~ l/a . The spread- 
ing during the recollision time t' rec can be estimated as 



Ax' ~ Ay' ~ Az' 



Ap' 



m 



(56) 



D. Muon production from an e + e plasma 

So far, we have studied the laser-induced process Ps — ► 
As was mentioned in the introduction, the cor- 
responding reaction rate we denote as the coherent rate. 
It is given in approximate form by Eq. (54). When N Ps 
atoms interact with the laser field, the total coherent rate 
R^P resulting from this ensemble increases proportion- 
ally with N, provided that each Ps atom independently 
creates a muon pair. The latter condition is fulfilled, 
when electrons and positrons originating from different 
Ps atoms do not interfere. This means that the electron 
wave packets from different atoms do not overlap, which 
is guaranteed according to the above, when a£\ <C n -1 / 3 , 
with the Ps density n. 

We have seen in Sec. III.C that the coherent rate 
is reduced by quantum mechanical wave-packet spread- 
ing of the initially well-localized electron-positron state. 
Spreading will be absent if, instead of forming a bound 
state, the electron and positron are initially free particles. 
In practice this could be realized by an e + e~ plasma. 
When driven by a strong laser field, e + e~ collisions occur 
inside the plasma, which can give rise to produc- 
tion as before. The corresponding reaction rate we call 
the incoherent rate in order to distinguish it from the co- 
herent rate resulting from Ps atoms. The incoherent rate 
scales with N + N-, where N + and 7V_ are the numbers of 
electrons and positrons, respectively. In a neutral plasma 
(7V + = iV_ = N), the total incoherent rate R^_ raises 
like N 2 . 

The incoherent rate is calculated via the square of the 
amplitude (3). Integrals like the one in Eq. (19) are ab- 
sent then. Along the same lines followed in Sec. III.C, 
one can derive an approximate expression for the incoher- 
ent rate of muon production from a nonrelativistic e + e~ 
plasma. It reads 



R P 



_L cf_ 

2 3 tt 2 m 2 <e 4 



2 

min 



N + N- 
V 



(58) 
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with the interaction volume V. In our situation, the lat- 
ter is determined by the laser focal spot size. 

Similar as before, Eq. (58) can be equipped with an 
intuitive meaning. To this end, we introduce the number 
of produced muons during the interaction time T and 
rewrite Eq. (58) as 



R,, 



V 



VT 



r^> (J- 



(59) 



Here, a is the field-free cross section of Eq. (55) and n± = 
N±/V are the electron and positron densities. The rep- 
resentation (59) has the advantage that the combination 
Nfj,/ (VT) is Lorentz-invariant and can directly be trans- 
formed into the cm. frame. The latter moves with the re- 
duced velocity (3 of Eq. (48) and the Lorentz factor 7 ~ £. 
The particle densities n'± in the cm. frame are related 
to the lab frame densities by n± = "f(n'± — (3j' z ) = ^fn' ± 
since the particle current density j' z vanishes. We can 
therefore recast Eq. (59) into the form 



Re+t 



Rl, 



v 



V 



or R' e+e - ~ oN + N- /V' . This means that in the cm. 
frame, the number of events per volume and time is re- 
lated to the corresponding cross section in the usual way. 

We can compare the incoherent rate with the coher- 
ent rate for a single Ps atom by setting N + = N- = 1 
in Eq. (58). The interaction volume in the lab frame is 
determined by the wave-packet volume V = (a£A) 3 /£ 
(see Sect. III.C), where a factor of l/£ accounts for the 
Lorentz contraction along the z axis. With these choices, 
the incoherent rate gives the same order of magnitude as 
the coherent rate in Eq. (57). 

To conclude this section we note that the coherent 
channel of muon production from Ps atoms evolves into 
the incoherent channel in the limit of large wave packet 
size (a£A > n -1 / 3 ). On the one hand, with increasing 
wave packet size the coherent rate is more and more sup- 
pressed; on the other hand, electrons (positrons) stem- 
ming from different Ps atoms start to overlap and the 
gas of Ps atoms transforms into an e+e~ plasma. There- 
fore, in the limit of large wave-packet size, the incoherent 
reaction rate (58) will eventually dominate the muon pro- 
duction process. 



E. The influence of the laser polarization 



In a previous paper, we have treated the process 
Ps — > in a laser field of circular polarization [12]. 

The threshold intensity in this case is larger by a factor 
of two, since the same minimum peak field strength is 
needed. From the investigation of laser-driven recollision 
processes in atoms or molecules (e.g., high- harmonic gen- 
eration) it is known that the efficiency of these processes 
is strongly suppressed in a circularly polarized laser field 
[45-47]. The reason is that the ionized electron is forced 



onto a circular trajectory and does not return to the par- 
ent ion, which stays at rest due to its heavy mass. In 
our case, the situation is different as both binding part- 
ners move along circular orbits (within the polarization 
plane) under the influence of a circularly polarized field. 
But still, we also find a strong reduction of the reaction 
Ps — > n + when the laser polarization is circular. The 
corresponding rate approximately amounts to 



5(c) _ <r { a 
Ps ~ al Uc 



4 / 4m \ 2/3 



(60) 



with £ c = ea/m. Here, the Volkov state normalization of 
Eq. (4) and the approximate treatment of the (5-function 
in Eq. (51) has been applied. We emphasize that this 
leads to a correction factor of £~ 2 as compared to Eq. (52) 
in [12], resulting in an even smaller production probabil- 
ity than predicted there. The physical conclusions of 
this paper thus remain unchanged. In particular, the 
rate suppression in a circularly polarized field is mainly 
expressed by the damping factor [ao/(A£ c )] 4 <~ 10~ 25 in 
Eq. (60). This factor can be attributed to the classical 
motion of the electron and positron, which corotate on 
opposite sides of a circle of radius <~ A£ c in the polariza- 
tion plane. The interparticle distance is thus larger than 
their wave-packet extensions. This is in contrast to the 
case of a linearly polarized field, where the classical tra- 
jectories periodically meet and the rate damping results 
from quantum mechanical dispersion (cf. Sec. III.C). 
Equation (60) can also be written in the form 



It 



(c) 
Ps 



15 a 2 



16 m 2 e c &a(&A) 3 \m?; 



1/3 



(61) 



which is more similar to Eq. (54) and allows an interpre- 
tation along the lines of Eqs. (55)- (57). According to the 
above, the electron-positron distance in the cm. frame 
amounts to Ax' ~ Ay' ~ A£ c , Az' ~ a£ c \ [cp. Eq. (56)]. 
Hence, as before in Eq. (57), the total reaction rate can 

be expressed as i?pj) <~ aj (£, c Ax' Ay' Az'), where we have 
dropped factors which are of order unity at £ c > M/m. 

The incoherent rate of muon production from a non- 
relativistic e+e~ plasma driven by a circularly polarized 
laser field reads 

R {c l 

ft' ft 



9 a 2 N+N- 

lOTT TO 2 £f V 



(62) 



The formula involves an additional factor of ~ £,r 2 in 
comparison with Eq. (59) in [12], as well. Since quan- 
tum wave-packet spreading and classical motion are not 
crucial here, the incoherent rate is of the same order of 
magnitude as the corresponding one in a laser field of lin- 
ear polarization [cf. Eq. (58)]. The incoherent muon yield 
in a circularly polarized field is slightly larger, which can 
be attributed to the fact that the field amplitude is con- 
stant and not oscillating, as it is for a lineraly polarized 
wave. Similar results are known from other strong-field 
phenomena such as e + e~ pair creation in combined laser 
and Coulomb fields (see, e.g., [48]). 
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IV. NUMERICAL RESULTS AND DISCUSSION 



A. Total muon yields 



On the basis of Eq. (22) we have also performed numer- 
ical calculations of the laser-driven process Ps — > n + ^~ ■ 
The high-order generalized Bessel functions were evalu- 
ated with the asymptotic expansion (30) . Figure 2 shows 
the dependence of the total reaction rate on the laser 
peak field strength for a near-infrared laser frequency of 
u) = 1 eV and a single Ps atom. For comparison, the fig- 
ure also contains the analytical rate estimate in Eq. (54) 
and the corresponding results for a circularly polarized 
laser field. For both polarization states, the agreement 
between the numerical simulation and the approximate 
expression is reasonably good. The reaction rates in the 
case of circular polarization are about eight orders of 
magnitude smaller than in the case of linear polarization. 
The very pronounced impact of the laser polarization on 
the process efficiency was discussed in Sec. III.E. 

In a linearly polarized driving field, the muon produc- 
tion rate starts from zero at the energetic threshold [see 
Eq. (54)] and reaches maximum values of the order of 
10~ 10 s _1 at laser intensities of about 5 x 10 22 W/cm 2 
(£ s=s 170). The overall behaviour of the rate is simi- 
lar to that of the cross section for the field-free process 
e+e~ — > n + in Eq. (55). At higher intensities, how- 
ever, the rate decreases more strongly than the cross sec- 
tion, since the decrease results not only from the typi- 
cal high-energy dependence 1 / E 2 m but also from the en- 
hanced wave-packet spreading [cf. Eq. (57)]. 
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FIG. 2: Total rates for the process Ps — > fj, [i~ induced by 
an intense near-infrared laser field (ui — 1 eV) , as a function 
of the laser peak field strength. The black squares and the 
solid line refer to a linearly polarized laser field and show the 
results of numerical calculations based on Eq. (22) and the 
analytical estimate in Eq. (54) , respectively. The black circles 
and the dashed line show the corresponding results for a laser 
field of circular polarization. 



For comparison we note that the rate for muon cre- 
ation in a hypothetical large-scale e + e~ collider experi- 
ment, that employs a single electron and a single positron 



only, would be of the order of 10 _20 s _1 [12]. The rate 
resulting from the "e + e~ micro-collider" based on Ps is 
orders of magnitude larger, due to the small impact pa- 
rameters and the correspondingly large current densities 
("luminosities") achievable with linearly polarized driv- 
ing fields. A real collider experiment, however, employs 
beams of ~ 10 10 particles and reaches reaction rates of 
about 1 s _1 this way. This number, in turn, is orders of 
magnitude larger than the coherent rate from a single Ps 
atom. In Sec. IV. C below, we discuss the implications 
of these rate estimates for an experimental realization of 
the laser-induced process Ps — > n + ^~ ■ 



B. Muon energy distribution 

In Fig. 3, partial production rates with respect to the 
number r of absorbed laser photons are shown (i.e., the 
contributions to the total rate stemming from a net ab- 
sorption of r laser photons in the production process). 
The photon number is given in units of tq = m£ 2 /(2o;). 
Since ruj rts P° [see Eq. (42)], the partial rates reflect 
the energy distribution of the created muons. The dis- 
tribution is broad, indicating a substantial width of the 
energy which is available in the process and transformed 
into muons. 

First we emphasize that the broadening of the muon 
energy spectra is not due to the initial momentum spread 
in the Ps atom, which amounts to Ap ~ ma. In ac- 
cordance with Eq. (55), the cm. energy associated with 
the e + e~ annihilation is E cra = [(q + + q_ — nk) 2 ] 1 / 2 ~ 
200 McV. It depends on the initial momentum p via the 
effective particle momenta q± [see Eq. (5)]. The ini- 
tial momentum spread Ap translates into a longitudi- 
nal width A(q+ + q-) z ~ ma 2 £ 2 , whereas the transver- 
sal momentum components add up to zero. Since Ap is 
small, the number n of photons absorbed at the electron- 
positron vertex of the Feynman graph varies within nar- 
row bounds [see Eq. (14)]. As a result, the width of the 
collision energy AE cm ~ E cm a 2 <~ lOkeV is small. 

The main limitation for a well-defined effective collision 
energy comes from the fact that the total number r of ab- 
sorbed laser photons is not fixed, but has a large width 
Ar (see Fig. 3). This is because the number of photons 
absorbed at the muon-antimuon vertex varies within a 
wide range. The contribution of these laser photons leads 
to an enhancement of the available energy, that can be 
transformed into muons. The effective collision energy is 
accordingly given by E c . m . = [(q + + g_ +rfc) 2 ] 1 / 2 > E cnl , 
as indicated by the S function in Eq. (21). The avail- 
able energy must be larger than E cm w 2M, because 
the muons have to be produced with their effective mass 
M* > M in the laser field. The energetic difference 
AM 2 = 2(M 2 — M 2 ) is supplied by the minimal num- 
ber of laser photons r m j n , according to Eq. (33). At 
the laser intensity threshold (where the process prob- 
ability vanishes), only photon numbers r w r mm con- 
tribute and the muons are produced with a sharp energy 
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P w M 2 /m w 22GeV. We have seen in Sec. III.B that 
in this case the muons are created in the field at rest in 
the cm. frame. For larger intensities, the range of con- 
tributing photon numbers quickly increases and attains 
a typical size of Ar ~ 10 10 for the parameters consid- 
ered in Fig. 3. This leads to a corresponding spread 
Aruo <~ lOGeV of the muon energies in the laboratory 
frame (or Arui' ~ 100 MeV in the cm. frame). The 
longitudinal and transversal momentum widths amount 
to AP Z ~ Arui ~ lOGeV and AP ± - Arw' ~ 100 MeV, 
where the first estimate results from the relation P z wF , 
and the second from U s=s r, in accordance with the prop- 
erties of the Bessel functions. 

The shape of the curves in Fig. 3 can be understood 
in more detail via the kincmatical analysis in Sec III. A. 
First, according to Eq. (33), the minimal number of laser 
photons required from kinematical constraints amounts 
to r m ; n = 2.2 x 10 10 , independent of the laser intensity. 
If we express this number with respect to the respec- 
tive values of r , we find r min /r = 3.0 (£ = 170), 2.2 
(£ = 200), and 1.4 (£ = 250), as is displayed in Fig. 3. 
The partial reaction rate for r = r m ; n is always zero. 
Furthermore, the curves exhibit maxima at r w f, in 
agreement with Eq. (36); the typical photon number f is 
increasing with the applied laser intensity. We point out 
that the width of the curves is significantly larger than 
in the case of circular laser polarization (see Fig. 3 and 
Ref. [12]). This can be understood by observing that, in 
a linearly polarized field, e + e~ collisions not only occur 
at the optimum laser phase (which is predicted by the 
simple man's model and gives rise to the absorption of f 
photons), but also at earlier or later phases. Hence the 
collision energy varies, which broadens the energy dis- 
tribution of the created particles. Different from that, 
in the circular polarization case all creation phases are 
equivalent, so that the width of the partial production 
rate is purely of quantum mechanical origin. 



C. Experimental feasibility 

The maximum coherent rate resulting from a single 
laser-driven Ps atom is of the order of Rp s <~ 10 _10 s _1 
(see Fig. 2). The total process probability is obtained by 
multiplying the rate with the interaction time, which is 
given by the laser pulse duration. The high intensities 
required to ignite the process (I ~ 10 22 -10 23 W/cm 2 ) 
are attained in short laser pulses, with typical durations 
in the ps-fs domain. The muon yield per laser pulse from 
a single Ps atom is therefore very small (< 10~ 20 ). 

Since the total probability of the coherent process in- 
creases linearly with the number of Ps atoms involved, 
high Ps densities are desirable. In recent years, a re- 
markable progress has been made in the efficient pro- 
duction, accumulation and trapping of Ps atoms [49]. 
The highest Ps density achieved so far is of the order 
of n = 10 15 cm -3 [50]. At this density, a typical laser fo- 
cal volume of V = (10A) 3 ~ 10 





l,2xl0" 21 - 


0> 


l.OxlO" 21 - 


OS 

c 


8,0xl0" 22 - 


luetic 


6,0xl0" 22 - 


proc 


4,0xl0" 22 - 


partial 


2,0xl0" 22 - 


0,0- 



I / 


\x2 


I / 






; s y x5 - 


■ x 10 s 




— ' ■ ' i — U 1 — 


i— * ■ 



5 10 15 

photon number r/r 



20 



FIG. 3: Partial rates for the laser-driven decay Ps — > /v n~ , 
as a function of the number of absorbed laser photons r. The 
linearly polarized laser field has a frequency of w = 1 eV and 
an intensity parameter of £ = 170 (solid line), 200 (dashed 
line, enhanced by a factor of 2), and 250 (dash-dotted line, 
enhanced by a factor of 5), respectively. For comparison, the 
dotted line shows the corresponding result for a circularly 
polarized field with £ c = 250; it is enhanced by a factor of 
10 8 . The photon number is given in units of ro = m£ 2 /(2ui). 



atoms. The total coherent reaction rate could accord- 
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ingly be increased by this number to R 
There are proposals how to generate even denser Ps sam- 
ples (n > 10 18 cm~ 3 ) with regard to the formation of a 
Ps Bose-Einstein condensate [50] or experimental studies 
on antihydrogen [51]. Such high densities, however, vi- 
olate the condition a£A <C n -1 / 3 for coherent collisions, 
so that electron (positron) wave packets from different 
Ps atoms will considerably overlap at the collision. The 
muon creation then proceeds via the incoherent chan- 
nel. Incoherent muon production from a laser-driven 
e + e~ plasma requires high target densities, as well, be- 
cause of the small focal spot size of strong laser pulses, 
which determines the interaction volume. By assuming 
again a typical focal volume of V = (10A) 3 w 10~ 9 cm 3 , 
which contains N± — 10 7 particles at a plasma density 
of 10 16 cm~ 3 [49, 51], we obtain a total incoherent rate 

of R ( e +l_ - 10~ 2 s- 1 at £ w 230 [sec Eq. (58)]. In view of 
the short duration of strong laser pulses and the above 
values of R^ and R^+l- > both the coherent and incoher- 
ent channels of laser-driven muon pair production seem 
to be hardly observable in experiment at first sight. 

The wave packet spreading, however, which has a detri- 
mental impact on the coherent reaction rate, can be 
strongly reduced by applying more complex driving-field 
geometries than a single laser wave. A suitable config- 
uration is formed by two counterpropagating laser laser 
waves, as is realized for example by the Astra Gemini sys- 
tem at the Rutherford Appleton Laboratory or the JETI 
"photon collider" at the University of Jena [52]. The re- 
sulting Ps dynamics has been investigated in Ref. [7] by 
means of a classical Monte-Carlo simulation for circularly 
polarized laser waves of equal frequency and intensity. It 
was found that the wave packet spreading is largely re- 
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duced since the recollision times are much shorter than 
in a single laser wave. Moreover, the resulting laser mag- 
netic field is oriented along the electric field and induces 
a focussing force [5] . As a consequence, the particle wave 
packets are as small as (10A) 3 at the recollision. This 
implies that the muon yield could be increased by <~ 10 
orders of magnitude as compared to the case of a sin- 
gle driving laser field [see Eq. (57)]. Moreover, due to 
the small wave packet size, higher Ps densities can be 
exploited, without violating the condition for coherent 
collisions. At the envisaged density of n > 10 18 cm~ 3 

[50, 51], a total coherent reaction rate of ~ 10 7 s _1 
could be achieved, with N — 10 9 . Under these circum- 
stances, the observation of coherent muon pair produc- 
tion from laser-driven Ps becomes feasible at high laser 
repetition frequency ^Hz-kHz [1, 2]. The resulting num- 
ber of events approaches ls^ 1 , which is comparable to 
the event rates at large-scale accelerators (see Sec. IV. A). 
As to the kinematics of the created particles, we note 
that the laboratory and cm. frames coincide here due 
to the symmetric setup. The expected emission pattern 
is therefore approximately isotropic, in contrast to the 
muon kinematics resulting in the case of a single driving 
laser field (see Sec. III. A). 

We point out that also the incoherent production 
channel is enhanced in the countcrpropagating beam- 
geometry since the cm. motion is absent. In accor- 
dance with Eq. (59), the enhancement factor amounts to 
£ 2 ~ 10 -10 5 . Another promising field combination is 
formed by a single laser wave and an additional static 
magnetic field, which is spatially confined [7]. This ge- 
ometry allows for higher collision energies, as the large 
longitudinal momentum (q z <~ m£ 2 ) of the electron and 
positron is transversally redirected by the magnetic kick 
and can thus be exploited for particle production. 

To conclude this section, we point out the following 
technical aspect. Our calculation assumes that the elec- 
tron and positron trajectories lie completely within the 
laser field. The transverse extent of the classical trajecto- 
ries along the laser electric-field axis is Ax <~ A£ ~ 10 2 A. 
This is larger than covered by a typical laser focal di- 
ameter. On the other hand, there is practically no mo- 
tion along the magnetic-field direction, apart from wave 
packet spreading (Ay ~ a£ A ~ A) . In order not to waste 
laser power, it would therefore be useful to employ laser 
foci of elliptical shape: prolonged along the polarization 
axis and compressed in perpendicular direction. This 
way, a laser focal area of order AxAy ~ (10A) 2 could 
fully accomodate the electron and positron motion in the 
field. 



V. SUMMARY AND CONCLUSION 

The creation of a muon pair by e + e~ annihilation in an 
intense laser wave of linear polarization was calculated. 
The initial electrons and positrons were assumed to form 
a gas of Ps atoms or a nonrelativistic plasma, whereas 



the energy required for the process is provided by the 
laser field. The minimum field intensity amounts to a few 
10 22 W/cm 2 in the near-infrared frequency domain, cor- 
responding to an intensity parameter of £ w 150. These 
laser parameters are within experimental reach [1, 2]. 

Our calculation proceeded within the framework of 
laser-dressed QED, employing the strong-field approxi- 
mation. The initial bound state was expanded into a su- 
perposition of products of Volkov states, weighted by the 
Compton profile of the Ps ground state. The produced 
muons were also described by relativistic Volkov states. 
The process amplitude has been evaluated in two ways. 
First, the standard method of expansion into Fourier se- 
ries was applied. The appearance of generalized Besscl 
functions of very high order ~ 10 11 represents a technical 
difficulty here. In a second approach, we evaluated the 
amplitude by the saddle-point integration method and 
obtained this way a suitable asymptotic expansion for 
the generalized Bessel functions, which has been used for 
numerical computations. The numerical results on the 
total reaction rates are in good agreement with corre- 
sponding analytical estimates. The latter allow for an 
intuitive understanding of the process in terms of the 
field-free muon production cross section and the size of 
the colliding electron and positron wave-packets. 

The produced muons are characterized by ultrarela- 
tivistic energies and very narrow emission angles around 
the laser propagation direction. The muonic life time is 
accordingly enhanced by relativistic time dilation. The 
peculiar muon kinematics is explicable in terms of a semi- 
classical simple man's model. It essentially arises from 
the fact, that the electron and positron are forced by the 
laser field to longitudinal velocities close to the speed of 
light, i.e., the cm. frame is moving at high velocity with 
respect to the laboratory frame. The energy distribution 
of the muons is broad, because the number of photons ab- 
sorbed from the laser field during the production process 
varies over a wide range. 

As to the total reaction rates, we found that a high 
scattering luminosity is achieved in the case with initial 
Ps atoms. This is due to coherent e+e~ collisions at 
small impact parameters, which result from the micro- 
scopic extension of the initial state and the symmetric 
electron-positron dynamics in the laser field. Limitations 
on the process efficiency are set by quantum mechanical 
wave-packet spreading which dilutes the particle density 
at the collision. The resulting coherent reaction rate from 
a dense Ps sample is about 10~ 4 s^ 1 for a laser of inten- 
sity w 5 x 10 22 W/cm 2 and focal volume w (10A) 3 . This 
value is several orders of magnitude larger than the cor- 
responding reaction rate in a circularly polarized field, 
where the process is suppressed by the classical electron 
and positron trajectories. The incoherent rate for muon 
production from a laser-driven e + e~ plasma is essentially 
insensitive to the field polarization and reaches values of 
~ 10~ 2 s _1 at high plasma densities. 

The field-induced process Ps — > is experimen- 

tally accessible in a setup that employs a high-density 
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2m 



Ps target (n <~ 10 18 cm" 3 ) and countcrpropagating laser where q±(0) denote the values of q± at p — and 
beams to control the wave-packet spreading. The co- 
herent production rate could be increased to ~ 10 7 s _1 

this way, which translates into an effective muon yield Aq(p) 
of about Is -1 when the laser pulse duration (~ 100 fs) 
and repetition frequency (~ Hz-kHz) are taken into ac- 
count. The event rate achievable with the laser-driven 
"e+e~ micro-collider" in the crossed-beams configuration 
is similar to those at large-scale accelerator facilities. An 
experimental observation of laser-induced muon pair cre- 
ation from Ps decay, though challenging, will therefore 
be rendered feasible by the next-generation of high-power 
laser devices. 



(A2) 



APPENDIX A 

In the derivation of Eq. (21) in Sec. II. A, we neglected 
weak dependences on the bound-state momentum p in 
the transition amplitude when performing the integral 
over the Ps momentum distribution. As regards kine- 
matic terms appearing as linear factors, this approxima- 
tion is certainly justified. But we also ignored the p de- 
pendence in the argument of the 5 function in Eq. (21), 
which is a more subtle simplification. In this appendix 
we show that this approximation is indeed justified. 

With the 5 function included, the integral in Eq. (20) 
reads 



d 3 p 

J^Kf 



x S(q^ 



$(p)J r (U-u,V-v) 
q--Q+-Q-+rk). (Al) 



According to Eq. (5), the effective momenta are weakly 
dependent on p (see also Sec. IV. B). We perform a Taylor 
expansion up to leading order 



<tf +<£*<#(0) + <£(0) + Ag(p) 



k^ 



Within the momentum range given by $(p), the correc- 
tion term can be as large as Aq(p) <~ ma 2 f/8. Since it 
is multiplied by the wave four- vector in Eq. (Al), we can 
consider Aq(p)/u) ~ 10 5 as a correction to the photon 
number r. Recall that a typical value is r ~ 10 11 so that 
r 1 / 3 <~ 10 4 . Hence, by the Bessel function properties, the 
correction to the photon number induced by the term 
(A2) could, in general, be important. 

By writing the square of the amplitude as a product of 
two terms, we can express the reaction rate as 



d?Q+ f d 3 Q_ 



d 3 p$(p)J r (U-u,V-v) 



x 5 (q+(0) + q_(0) + Aq{p)k/uj - Q + - Q_ + rk) 

x J d 3 p' §{p')J r {U -u',V- v') 

x S (q+(0) + <z_(0) + Aq{p')k/uo -Q+-Q-+ r'k) 

(A3) 

where linear kinematic factors were suppressed in the no- 
tation. Both 5 functions in Eq. (A3) contain the term 
<5(Q+ + Qf )• One of these we exploit to integrate by Q 1 ;, 
and the integral over Q + z is performed with the help of 
<*(9+,*(0) + g-,*(0) + Aq(p>) - Q +tZ - Q_, z + r'oo). The 
latter <5 function fixes the value of to 

Q+,z = q+, z (0) + q-, z (0) + Aq(p')-Q_ iZ +r'LO. 

Another two of the remaining S functions coincide and 
we obtain (with the interaction length L) 



L3 [ T^E fd 3 P $(p)Jr(U-u,V-v)6(q + (0) + q _(0) + A q (p)-Q + -Q°_+ru;) 

— r y 

x [ ijf$(p<)J r (U -u',V- v')6 (g» (0) + q°_(0) + Aq(p') - Q% - Q°_ + r'w) 

2 

L ' 3 / -7W- Yl f ~7== *(P) J r(U ~u,V- v)5 (q° + (0) + q°_(0) + Aq(p) Q% Q°_ + ru) 

J Q - r J JQ° + V ' 



(A4) 



where in the last step we made use of the fact that Q+. z 
and thus Q° + can also be expressed in terms of p and r, 
since Aq(p)+rui = Aq(p') + r'uj. By performing a Taylor 



expansion 



Q+ = JM* + Q\ + Q 2 +iZ w Q° + (0) + Aq(p) , 
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we observe that the leading p dependence drops out of 
the argument of the 6 function. Hence, we arrive at the 
desired result 

Rps ~ L i q^o) ^I/ <pp*w j *( u - u > v - v tf 

x 6 (q° + (0) + g° (0) -Q%- Q°_ + rw) , 

which shows that the argument of the 5 function in 
Eq. (21) may be taken at p = 0. 



APPENDIX B 

In this appendix we calculate the integral in Eq. (19): 
d 3 P z 



Jr = J Hp)MU -u,V-v), (Bl) 



where 



„-£( 1 + d-i) (B2) 

according to Eq. (11). In the Taylor expansions on the 
right-hand side of these equations we kept those terms 
which are of importance according to the properties of 
the function J r (U — u, V — v). By virtue of Eq. (8) we 
can expand Eq. (Bl) into series J r — J2 S Jr,s with 



Jr.s 



d 3 P - 



(2tt) 3 



$(p)J r -2s(U -U)J S {V -V). (B3) 



From the kinematical analysis in Sec. III. A we know that 

V w v. As a consequence, the magnitude of the argument 

V — v is much smaller than the order r of the generalized 
Bessel function. This is corroborated by our numerical 
calculations, according to which \V — v\ is smaller than 
r by an order of magnitude in the integration region of 
interest. In order to obtain significant values of J r ^ s it 
is required that s w V — v and r — 2s w U — u. Since 
s <C r and u ~ £,p x /u ^ ^am/uj <C 17 ~ £ 2 m/oj, this 
implies r w £/ '. Now, J r -2s(U — u) is almost constant 
for |u| < Umax ~ and strongly damped otherwise. 
Therefore, we can perform the p x integration in Eq. (B3) 
approximately and obtain 

J r , s « 2p£ u »J r _ 2i ([0 / %% *(p) J S (V - v) , (B4) 



(2tt) 3 



where 



P" 



io ( m£ z 



2\ 1/3 



is the p x value corresponding to u max by Eq. (B2), and 
<&(p) may be taken at p x — since p™ ax <C ma. By 



introducing the variables a = p y + p z and r = p 2 — p 2 , 



Eq. (B4) becomes 

2p ma ; 



Jr.s 

with 



(2^ 



poo 

J r -2s{U) da$(a)I s (a) , 
Jo 



-a \fo 

L 



: J, 



V-vo(l 



i 1 2m 2 ) 



AV+C 



J s (x)da 



(B5) 



Iav-c VC 2 - {x - At/) 2 ' 

where the substitution x = V — v [l — r/(2m 2 )] was made 
in the second step and the definitions ( = —avo/(2m 2 ), 
AV = V — vq were employed. In Eq. (B5) the width 
of the function J 8 {x) is determined by s 1//3 ~ AV 1 / 3 , 
whereas the width of the square root in the denominator 
is given by £. In order to evaluate the integral I s (a), we 
accordingly divide the integration range into regions and 
show 



( 7rj s (Av) (c«Ay 1 / 3 ), 

I (C>AV 1/3 and| S -AV|«C), 



(C> AV 1/3 and \s- AV\ ~ £). 



(B6) 



In the first range (( < AV 1/3 ), we have J s (x) w J S (AV) 
and thus obtain 



« J S (AV) / 



Ay+c 



av- ( VC 2 ~(x- AV) 2 
= nJ s (AV). 

The second range (C > AV 1/3 , \s - AV\ < C) similarly 
yields 

7i 2) M« ? / J s (^x«-. 

In the third range we employ the asymptotic expansion 
[42] 



1 2s/2 (3-x) 3 / 2 



and find 



i r s 



dx 



K 
r(|) 



2V^ (s-x) 3 ' 
3 



3 2/ 3 y^T75 ^CAVVa ' 



where T denotes the Gamma function. This shows 
Eq. (B6). 
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Now we can evaluate the integral I s 
da $(<r)I s ((j) over a in the corresponding regions: 



2m 2 

I «o| Jo 
2m 2 



Pol 



dC$(<7)7rJ s (AV) 
AV 1/3 $(0)7rJ s (AV), 



(B7) 



where $(cr) was approximated by $(0), since the range 
of £ values corresponds to cr <C m 2 a 2 . With the abbre- 
viation y = 2m 2 aQ/|w |, we find in the second region 



li 2) 



2m 2 f°° d( $(0) 

bo I 7 AV V3 C (i + yC) 2 

2m 2 



Pol 



$(0) In (l + yAV 1/3 ) - In (yAV 1/3 ) 
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Pol 
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l + 3/AyV3 
/ a 2 \v \ 
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(B8) 



where yAV 1 / 3 = 8AF 1 / 3 /(a 2 |w |) < 1 was exploited in 
the last step. Note that we have not restricted the range 
of integration according to \s — AV| <C C here, since 
1/C *C (AV 1 / 3 . For calculating the third region, 
we do restrict the integration to satisfy \s — AV| <~ ( as 
follows: 



2m 2 



2 r \s-AV 1 ' 3 \+AV 1 / 3 



$(CR 



Pol J| s -Ayi/3|_ A yi/3 y/^AV 1 / 3 



- AV| - AVV3 



2*(0)i 



AV 1 / 3 



AV| 



(B9) 



where in the last step AV 1 / 3 <C C ~ l s ~ AV| was applied. 

according to Eqs. (B7)-(B9), 



Since 7 S (3) < /.. 



r(2) 



r(i) 



we may write i s w 2j| 1 - ) and find 
m 2 p 



-,2^max /\"\/l/3 



■$(0)- 



|t) | 

By summing over s, we finally obtain 



Jr-2s{U)J s {AV). 



$(0). 



m lo 



2 s m e\ 1/3 av 1 / 3 



Pol 



2v / 2tt 2 £ 
which coincides with Eq. (20). 
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APPENDIX C 



In this appendix we derive the asymptotic expansion (30) 
for the generalized Bessel function J^(U,V) (see also 
Refs. [31-33]). We consider the integral J = Jdy-e if( - y ~^ 
of Eq. (25). Due to periodicity and the unbounded range 
of integration, Eq. (28) defines an infinite number of 
saddle points, which can be divided into four classes 
(j = 1, 2, 3, 4) according to 



(CI) 



Here, N is an integer number and the four particular 
solutions 



toy 
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(2) 



anros ( + p ] . 

( u 

2tt — arccos — — — + p 
V 8 V 



ivy^ = arccos ( -— - p 



ujy^ — 2ir — arccos ( — ^ — p 



(C2) 



with 



U 2 AQ 1 



64V 2 Wuj 2 ' 

lie in the intervall [0, 2tt). At the saddle points, the sec- 
ond derivative of the phase (26) becomes 



dy 2 _ 



= ±8u 2 Vpsm{ujy { !' N) ) 



(C3) 



with the positive (negative) sign for j = 1,2 (j = 3,4). 
Using the general identity 



E « 

N= — oo 



-2,rijV=g 



w £ c5(AQ-Au), 



iV=— oo 



we can sum up all contributions from the saddle points 
y^' N ^ in closed form and obtain 

oo 

T« = V^ 1 S(AQ-Nlj) 

iV=— oo 

iNiuy W -iU sin(ujy W )-iV ein(2wy !_ 1) )+CT< 1 )if 



^/8p|Vsin(w2/i 1 



where cr' 1 ) = sgn[V sm(ujy^)] fixes the sign of the second 
derivative in Eq. (C3). Analogously, we find 1^ = jW* 
so that 

oo 

j(i)+j(2) =2V2^ ^ 6(AQ-Nu) 

N=—oo 

p iArwy ( i 1) -it/sin(wiy ( i 1) )-iysin(2wy ( i 1) )+(T (1) i^ 

xRe 
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Note that sm(wy^) = - sin(wyi^). The contributions and we finally obtain 
of the saddle points y^' N ^ and y^' N ^ add up similarly 



J = T« + 1^ + 1^ + T< 4 ) = 2V2^ ]T S(AQ - Nlo) 



x Re 



N=-oo 

e iAfwy ( L 1) -i(7sin(wy ( L 1) )-iysin(2wj/ ( L :L) )+ ( T <1) if e i7Vwiy ( L 3) -i;7sm(wy ( L 3) )-iVsin(2w 2 /^ ) )-<j (3) i 



+ 



(C4) 



with o^ 3 ) = sgn[V r sin(wy_')]. A comparison of Eq. (C4) 
with the last line of Eq. (29) confirms the asymptotic rep- 
resentation (30) for the generalized Bessel function. 

We comment on the applicability range of Eq. (30) [32]. 
First of all, in the spirit of the saddle-point method, large 
values \U\, \V\, \N\ ^> 1 are required. The saddle points 
(C2) should be real and well separated along the real axis. 
This implies p 2 > and | cos(u;y_)| < 1. If not all saddle 
points are real, then the contribution from the complex 
ones may be ignored to a good approximation, as illus- 
trated below (cf. also Eq. (3.7) in [32]). Moreover, the 
value of the second derivative d 2 f/dy 2 _ in Eq. (C3) must 
not be too close to zero; otherwise the next-order term in 
the Taylor expansion of the exponent (26) becomes dom- 
inant. A comparison of the asymptotic formula (30) with 
the exact evaluation of Jn{U, V) via the series expansion 
(8) is displayed in Fig. 4. We emphasize that only real 
saddle points have been taken into account here; contri- 
butions from saddle points with non-zero imaginary part 
are exponentially damped. I.e., in this case only one of 
the two terms within the real-part function in Eq. (30) 
was considered. We see that the asymptotic expansion 
provides a very good approximation, in general. Only 
in a small region around U = 15, V = —15 in Fig. 4(a) 
[or U = 80, V = —15 in Fig. 4(c)] the asymptotic for- 
mula fails. As shown in Fig. 4(b), it is exactly at this 
point where d 2 f/dy 2 _ vanishes. In order for Eq. (30) to 
be applicable, the parameters need to fulfill the condition 

Here, /" = ±8Vpsm(uy_) and /'" = 8F[±pcos(o;j/_) - 
sin 2 (wy_)] denote the second and third derivatives of / 
with respect to ojy- , taken at the potential saddle-point 
positions u>y^ according to Eq. (C2). When /" tends 
to zero, the condition (C5) will be violated eventually, 
since /"' remains finite [see Fig. 4(b)]. In this situation, 



moreover, | cos(wj/_)| tends to unity, so that the two re- 
spective saddle points approach each other. Then the 
isolation criterion is no longer satisfied. 

In our numerical calculations the asymptotic expan- 
sion (30) is applied to evaluate J r (U, AV) in Eq. (20). We 
check that the above conditions of applicability are sat- 
isfied in the integration region of interest. In accordance 
with the kinematical analysis in Sec. III. A, the charac- 
teristic magnitudes of the parameters are r w \U\ <~ 
m^/uj ~ 10 10 and \AV\ - |V|/10 w m^/(20u) - 10 9 . 
The parameter p is therefore positive and can be approx- 
imated by a Taylor expansion as 

\U\ f 8rAV 16AV 2 \ 

9 ~ 8|AV| \ + U 2 + U 2 J 

\U\ , r 2\AV\ 
~ 8\AV\ \U\ + \U\ ' 

where the first term typically is larger than unity, the 
second term is close to unity, and the last term is much 
smaller than unity. According to Eq. (28), there are real 
saddle points z located at 

.c^™^-^, (C6 ) 

From a physical point of view we note that a cosine value 
close to unity agrees with the discussion in Sec. III.B, 
since the electron-positron annihilation occurs close to 
the field maxima (i.e., at | cos(u;a;_)| w 1) and the 
muons are created shortly afterwards. From Eq. (C6), 
we obtain |sin(u>z)| w 2^\AV/U\ < 1. This implies 
| /"| ~ | /'"I - 8\AV\ > 1 and the condition (C5) is 
satisfied. As a result, the asymptotic expansion (30) is 
applicable for our purposes of numerical computation in 
the mainly contributing regions of the integration domain 
in Eq. (22). 
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FIG. 4: Illustration of the applicability range of the asymp- 
totic expansion in Eq. (30) for the generalized Bessel function. 

(a) J$o(U, 15) in the range < U < 90; the thick dashed line 
shows the exact value resulting from the series expansion in 
Eq. (8) , the thin dotted line shows the asymptotic formula. 

(b) The condition (C5) for the example of the figure in (a), 
with the thick dashed and thin dotted lines showing the left- 
and right-hand sides of Eq. (C5), respectively, (c) Same as 
(a), but for J 50 (80, V) in the range -40 < V < 20. 
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